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A polymer chain pinned in space exerts a fluctuating force on the pin point in thermal equilibrium. 
The average of such fluctuating force is well understood from statistical mechanics as an entropic 
force, but little is known about the underlying force distribution. Here, we introduce two phase space 
sampling methods that can produce the equilibrium distribution of instantaneous forces exerted by 
a terminally pinned polymer. In these methods, both the positions and momenta of mass points 
representing a freely jointed chain are perturbed in accordance with the spatial constraints and the 
Boltzmann distribution of total energy. The constraint force for each conformation and momentum 
is calculated using Lagrangian dynamics. Using terminally pinned chains in space and on a surface, 
we show that the force distribution is highly asymmetric with both tensile and compressive forces. 
Most importantly, the mean of the distribution, which is equal to the entropic force, is not the 
most probable force even for long chains. Our work provides insights into the mechanistic origin of 
entropic forces, and an efficient computational tool for unbiased sampling of the phase space of a 
constrained system. 


I. INTRODUCTION 


According to the second law of thermodynamics, a system tends to seek a higher entropy state. When the increase 
in entropy can be gained through a change in a spatial coordinate, an effective force emerges. This force, purely 
entropic in origin, is a universal phenomenon, not specific to the underlying microscopic Hamiltonian of the system. 
Examples of entropic force are found in all realms of physics such as the Casimir force [T] , elastic force, and depletion 
force [5] . The concept of entropic force has even been applied to gravity mi]. 

Polymers provide an excellent model system to study the entropic force both theoretically and experimentally. They 
can be attached to a surface, confined in a volume, or constrained in a particular shape. Compared to polymers free 
in solution, geometrically confined polymers have reduced entropy. This reduction in entropy leads to an entropic 
force, which may influence the rate of detachment from a surface or the escape time of a polymer from confinement. 
The confinement-mediated force was measured for a chain in a tube [3, and calculated using semianalytical means 
and brownian dynamics simulations [B]. A simple heuristic argument estimates its magnitude to be ~ ksT/a, where 
a is the monomer length. 

More involved calculations include computation of the full partition function under confinement and its derivative 
with respect to displacement [7]. Alternatively, this force can be calculated from the local concentration of mass 
points near the surface mill- For a chain tethered to a hard wall, the force increases with chain length, but saturates 
at ~ ksT/a. For a temperature 3I0K and a monomer length of 3.4 A, this corresponds to 12.6pN. This is consistent 
with 13.3 pN value for the pressure integrated over the entire surface found by Bickel et al. BH], where a shorter 
monomer length of 3 A is assumed. 

These methods are adequate to calculate the entropic force as the mean thermodynamic force conjugate to displace¬ 
ment. However, entropic force itself is induced by a fluctuating quantity m- Several examples in biology point to 
the importance of rare events that deviate from the average behavior mis]. Likewise, force fluctuations may have a 
significant impact on chemical and biological processes [Mj. Therefore, finding the full distribution of instantaneous 
forces or its higher moments can give insights not easily accessible by statistical mechanics. 

In this study, we obtain the equilibrium force distribution due to a flexible polymer pinned in space or to a surface 
(Fig.[T]). The key idea of our approach is to make random moves in the phase space of generalized position coordinates 
and their conjugate momenta and calculate microscopic forces exerted by the system. This approach is different from 
conventional Markov chain Monte Carlo methods that typically apply moves only in the position space. Because of 
the universality of entropic force, we expect our method to be applicable to a wide range of dynamic systems. 
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II. METHODS 


A. Entropic force prediction 


The entropic force can be found from a thermodynamic calculation 


F = 


dA 

dr 


-keT 


dlogQ 

dr 


( 1 ) 


where A is the free energy, Q is the partition function, and r is the generalized coordinate of interest. By measuring 
the change in the partition function as the chain is slightly displaced from the pin point, we can compute the entropic 
force. As shown in Fig. Hb), when the tether is lengthened, additional conformations become available. For the full 
space case, the number of additional conformations is simply proportional to the area available to the first mass point. 
For the half-space case, we count the fraction of chain conformations compatible with the boundary condition m, 
and calculate the small increase (AQ) upon varying ri by Ari to numerically approximate Eq. 


B. Phase space sampling 

1. Monte Carlo sampling 

Beginning from some initial conformation and initial generalized momenta, we can generate a set of conformations 
for a coarse-grained chain in a heat bath. The Hamiltonian of the chain in terms of generalized positions (q) and 
momenta (p) will be given by 


i?(q,p) = C/(q) + ^P^M-ip (2) 

where the first term {7 is a potential energy based on chain deformation, and the second term is the kinetic energy 
expressed with the mass metric M. For a freely jointed chain (FJC), U is zero. We use the FJC model in this study 
because we want to obtain forces purely entropic of origin, but not due to physical interactions. The FJC model 
effectively coarse-grains an arbitrary polymer into Kuhn-length segments. The mass metric can be found from the 
Jacobian 

(M)y = (3) 

OQi 

M depends on q in non-Cartesian space. We choose 6 and (p angles of each link measured in the global reference frame 
as the generalized coordinates (Fig. [^c)). These are chosen over internal coordinates, as a change to the orientation 
of one link in the global frame will only affect two columns of the Jacobian matrix. 

At each step, one of the momentum coordinates can be perturbed, and the new kinetic energy computed. The trial 
step can be accepted or rejected according to the Metropolis criterion on the total energy. The position coordinates 
can also be perturbed, which will require computing the change in kinetic energy (via a change in M(q)) before the 
move can be accepted or rejected. 


2. Equilibrium sampling 

For longer chains, we sample points in the phase space using modal velocity decomposition m- The key idea is 
that microstates specified by both q and p can be directly sampled from the probability density function (p(q, p)) to 
build the distribution of any microscopic variable. Our method is different from methods categorized as equilibrium 
sampling m in that we sample both positions and momenta. For a canonical ensemble, /9(q, p) is given by 

p(q, p)d‘^^qd‘^^p cx (4) 

where j3 = l/fc^T. In this form, sampling individual variables and qi is not straightforward because they are 
coupled through M(q). Matrix diagonalization can be used to separate variables, but is not computationally efficient. 
Instead, the Cholesky decomposition can be used to express the symmetric, positive mass metric as 

M = mpfjJ 


( 5 ) 
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where /x is a lower triangular matrix. By defining modal velocities as v = /x^q, we can express the probability in a 
more separable form 

p(q, p)d^^qd‘^^p oc Vdet(M(q))d2^g • e" §(g) 

In this form, modal velocities (vi) can be sampled independently of each other from a Gaussian distribution. How¬ 
ever, positions (g^) become coupled through det(M(q)), which is known as the Fixman correction [TB]. Hence, we 
adopt a hybrid approach where we first sample a microstate using normally distributed modal velocities and weight 
the microstate-dependent variable by the Fixman term. Calculating the determinant of the dense metric tensor is 
computationally intensive, however Fixman demonstrated that this is equivalent to the determinant of a smaller, 
tri-diagonal metric of the constrained coordinates (for details, see Supplemental Information). 


C. Computing constraint force from positions and momenta 


We consider a flexible chain terminally pinned in space or to a surface. We can use Lagrangian dynamics to calculate 
the constraint force as a function of positions and velocities of the N mass points. Because distances between all 
adjacent mass points are constant, there are 2A^-|-1 degrees of freedom. We use three Cartesian coordinates {xi,yi,zi) 
for the first mass point, and 2{N — 1) angles {6 and cj)) for all other mass points of the chain (Fig. (c)). Choosing 
generalized coordinates this way simplifies subsequent calculations. The Lagrangian is given by 


L= -q'Mq 


( 7 ) 


For the first mass point of the chain to be radially constrained at distance r, we introduce a Lagrange multiplier A 


L' = ^q'^Mq -f A(ri - r) 


where ri = y/x^ +yl + zf. Using U in the Euler-Lagrange equation with respect to qk=i, 2,3 

d I dMi. . , dr I 


( 8 ) 


( 9 ) 


The first term on the right hand side vanishes as M does not depend on the Cartesian coordinates of the first mass 
point. Carrying out the differentiation on the left hand side produces an expression 


1\X •• I ■ ■ \dk 

Mkiqi + - 7 : - qiqj = A — 

oq, ri 


( 10 ) 


The second derivative of these terms must be solved for using the equation of motion for the angular coordinates. 

AX ■■ , • • 7 dMjk . . 

Solving for qj gives us 


.. f . . dMu , 1. . dMik\ 

Qi = I -qiqk— h -Qiqk— I 


dj 


dqk 


2 dqi J 


( 12 ) 


We can express this more concisely as the geodesic equation using the Christoffel connection coefficient, and making 
use of the symmetry in the first term under exchange of I and k 


qj = -qiqk^' 


Ik 


where 




dMki dM. 


V dqk 


-k 


jk 


dqj dqi 

Substituting this back in produces the Cartesian components of the constraint force 


qk . . f dMik 
A— = qiqj 
ri 


V dqj 

where Vi is the covariant derivative in the constrained subspace. 


— r\jMik^ = (jiqjViMjk 


(13) 


(14) 


(15) 
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III. RESULTS 


A. A single point constraint in space 


We first considerd a thermally equilibrated flexible chain with one end pinned to a single point in full three 
dimensional space. By using the FJC model without a potential term, we can focus solely on entropic contributions 
to the force. We used the Monte Carlo phase space sampling method to build a canonical ensemble of chains. As shown 
in Fig. [^a), kinetic energy is not equally partitioned between the mass points of the chain due to the constraints. 
The mass point nearest to the pin point has ~ 25% less kinetic energy whereas the endmost mass point has ~ 25% 
more. Most mass points in the middle have ksT on average, which would be expected for two degrees of freedom. 
Using these chain conformations and momenta, we calculated individual constraint forces from Lagrangian mechanics 
(Eq. |15[ ). The distribution of constraint forces is highly asymmetric, including both positive and negative values with 
a sharp cusp. The mean constraint force is a negative pulling force, and does not depend on the chain length. Due 
to high skewness of the distribution, the mean force 2kBT/a) is larger in magnitude than the most probable force 
IksT/a). The standard deviation also has little dependence on the chain length (Fig. [^d)). 

We compared this mean constraint force with the entropic force {F) predicted by statistical mechanics. To seek 
more conformations, i.e., to increase entropy, the chain would prefer detachment from the pin point. As a result, the 
chain exerts a radial entropic force on the pin point. If the distance (r) between the pin point and the first mass 
point increases by Ar, the number of conformations will increase in proportion to the spherical surface spanned by 
the mass point. Therefore, the partition function (Q) is proportional to 47rr^, and the force at monomer distance a 
is given by Eq. 

ksT Q{a + Ar) - Q{a) 2kBT 

-a;-= ~ 

This entropic force is identical to the mean constraint force in magnitude, and is also independent of chain length as 
shown in Fig. [^c). 

Interestingly, the mean constraint force due to the entire chain is equal to the mean centrifugal force exerted by a 
single mass point. A radially constrained mass point has only two degrees of freedom, and therefore, according to the 
equipartition theorem, its mean kinetic energy is equal to /csT. The mean of the centrifugal force (/c) is 

mv^ 
a 


2kBT 


(17) 



Since the energy {E = tou^/ 2) of a canonical ensemble is Boltzmann-distributed, the centrifugal force is exponentially 
distributed. 

p{E)dE = e-f^^dE = (18) 

As shown in Fig. [^b), this exponential distribution (dashed line) can explain much of the distribution of negative 
constraint forces. However, positive forces and a significant fraction of large negative forces deviate from this distri¬ 
bution. This result suggests that freely jointed mass points beyond the nearest one can have a collective impact on 
the pin point. 

To understand combined force generation by freely jointed mass points, we examined conformations that produce 
anomalous force values outside the range attributable to a single monomer. As shown in Fig. instances where the 
force has the positive sign result from sharp bending and large velocities. For the simple case of a trimer in two 
dimensions, we can derive a condition for compressive force 

elcos{e2-ei)>2el ( 19 ) 

This corresponds to the case where the centrifugal force exerted by on mi towards the pin point is more than 
twice that exerted by mi on the pin point in the opposite direction. 


B. A single point constraint on a plane 

We next considered a flexible chain attached to a surface. This confinement geometry arises in biology problems [HI 
I20j and polymer applications |21j . The phase space sampling method applies the same way except that we only accept 
chain conformations that do not cross the boundary. We computed constraint forces at different chain lengths. As 
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shown in Fig. |^a) and (d), the distribution of constraint forces is highly asymmetric, similar to the full space case. 
Most forces are exerted in the outward direction, but a significant fraction of forces are compressive. The mean 
constraint force increases slightly as a function of chain length up to ten monomer lengths, but quickly reaches a 
constant value. The mean constraint force plateaus at about 2.5kBT/a, which is 25 % higher than the full space case. 

For longer chains {N > 50), phase space sampling becomes computationally expensive. We therefore adopted an 
equilibrium sampling method to calculate constraint forces. This method requires weighting by the Fixman correction. 
We computed a raw force distribution (uncorrected, ‘x’ in Fig. [^, and later applied this correction (‘o’ in Fig. |^. 
Applying the Fixman correction lowers the average value of the force, as the largest forces are correlated with the 
largest extensions, which have a low weighting factor. This equilibrium method after correction produces similar 
distributions to the Monte Carlo method, highly asymmetric with more frequent negative than positive forces. 


To compute the corresponding entropic force, we consider the increase in the accepted chain conformations upon 
increasing the distance to the first mass point. The partition function can then be numerically differentiated using 
Eq. 16 to obtain the entropic force. The result did not depend on the size of displacement chosen. As shown by the 
red dashed line in Fig. |^b) and (e), the entropic forces based on conformational space only are slightly larger than 
the mean constraint forces {2.6kBT/a vs. 2.5kBT/a). 


IV. DISCUSSION 


We computed constraint forces exerted on a FJC in full space or half space by classical and statistical mechanical 
methods. Statistical mechanics was used to describe thermal fluctuations of the chain microstate, and classical 
mechanics was used to calculate the force required to constrain the chain against these thermal fluctuations. We 
found that the mean of constraint forces is equal to the entropic force calculated from the spatial derivative of the 
free energy. However, because constraint forces are asymmetrically distributed, the mean force is not equal to the 
most probable force. The force distribution is largely independent of the chain length beyond ten monomer lengths. 

In this study, we treat the FJC as a canonical ensemble, which exchanges heat with the environment at a constant 
temperature. Our equilibration of microstates in the phase space is similar to the Andersen thermostat |22| used in 
molecular dynamics with the added capability of simultaneously satisfying geometric constraints. At any instant, the 
mass points of the chain move with certain velocities at certain positions under the constraints set by the FJC model 
and the boundary. At a different time point, the mass points will adopt a different conformation with a different set 
of velocities as a result from collisions with the surrounding molecules of the heat bath. For each set of positions and 
momenta reset as a result of heat exchange, we calculate the microscopic constraint force. 

The most nontrivial result from this study is that the force distribution is not Gaussian even in the long chain limit. 
Because the mass points are constrained in the freely jointed chain model, their impact on the tether are neither 
identical nor independent. As a result, the constraint force does not obey the central limit theorem. Consequently, 
the most probable value and the average value of force do not coincide with each other. Therefore, the entropic force 
that we derive from statistical mechanics might not be the most relevant quantity in all circumstances. 

The fluctuation or second cumulant of an extensive thermodynamic variable can be easily calculated from the 
partition function when its conjugate intensive variable is fixed |2Jj . Textbook examples are energy fluctuation 
in a canonical ensemble and particle number fluctuation in a grand canonical ensemble. On the other hand, the 
fluctuation of an intensive variable such as force cannot be derived from the partition function when the conjugate 
extensive variable is fixed, i.e. in a fixed-distance ensemble [MllSi]- Here, we use the FJC terminally constrained in 
full space to illustrate this point. 

The constrained partition function (Qc) is given by (See Supplemental Information.) 


Q^(r) = J 


(20a) 





(20b) 


where He is the Hamiltonian of the constrained chain. The mass metric M is a 2N x 2N matrix that depends on 
2N global angles q. The dependence of the partition function on the distance between the end and the pin point 
(r) is explicity shown. Entropic force {F) arises because the constrained partition function increases with r, which is 
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reflected by the thermodynamic relation: 


F = 


1 dlogQcjr) 
P dr 

f 1 

J P dr 


d^^qd^^p/Q,ir) 


5i?c(q,p)e ,2N ,2N 


dr Qc(r) 

dHc(q,p )' 
dr 




(21a) 

(21b) 

(21c) 

(21d) 


The quantity in the constrained ensemble average can be identified with a microscopic force, which we denote as /. 
The second moment of / can be shown to be 




dr'^ 


/ d^H, 

\ dr'^ 


( 22 ) 


Whereas the first term can be calculated from the constrained partition function, the second term depends on the 
microstate of the system. If the rigid constraint is relaxed to a spring, the second term is simply the stiffness (/c) of the 
spring. In this case, force fluctuations can be easily calculated, similar to those for a particle trapped in a harmonic 
potential [MIEZ]- However, the force fluctuation diverges in the stiff limit (k —>■ oo) [28]. Therefore, to obtain force 
fluctuations with a hard constraint, one has to resort to computational means. 

We showed that for a terminally pinned chain in full space, our phase space sampling method produces constraint 
forces whose mean is equal to the entropic force 2kBT/a. Using Eq. 20b the entropic force can be worked out 


analytically. Integrating out the momentum dependence of Qc{f), one obtains 


Qcir) = const X J y/det(M)(i^^g 


(23) 


Qdr) depends on r through the mass metric M. It can be shown that •\/det(M) oc (Supplemental Information), 
and therefore the entropic force is given by 


F = - 


2kBT 


(24) 


In comparison, the entropic force due to the pinned chain in a half space is difficult to derive because the limits of 
integration imposed by the half plane depend on q in some complex fashion. 

People have used the method of images to calculate a simpler partition function in conformational space only (TJ [53] , 


Q{h) oc / d 

J half space 


2N 


q oc erf h\ 


2NaP 


(25) 


where h is the height of the pin point from the surface. By taking the derivative of \ogQ{h) with respect to h, a 
vertical entropic force can be derived. This vertical force increases with chain length for short chains, but is nearly 
constant at kBT/h for long chains. This force, however, is expected to be different from the mean constraint force 
obtained from our phase sampling methods for several reasons. First, the method of images yields only the vertical 
component of the entropic force out of mathematical convenience. Second, the analytical expression (Eq. 251 is correct 
only in the Gaussian limit of long chains. Third, the conformational partition function in Eq. 25 does not include the 
kinetic effect. However, it is unclear if removing the kinetic effect would significantly alter the mean force prediction. 
Therefore, we numerically computed the conformational partition function as a function of r instead of /i, and obtained 
a thermodynamic force conjugate to r using Eq. |16| This force is only slightly larger than the entropic force obtained 
by phase space sampling as shown in Fig. Bb) and (e). 

As the chain length becomes longer, the Monte Carlo method becomes computationally expensive. A faster method 
is to sample conformations directly from an a priori distribution. For example, to simulate a canonical ensemble of 
ideal gas, one can randomly sample velocities from the Maxwell-Boltzmann distribution and assign them to randomly 
positioned particles. Such equilibrium sampling, however, is not straightforward for a constrained system due to the 
coupling of p through M. By decomposing M, we can obtain modal velocities that individually obey a Gaussian 
distribution. Such a transformation skews the distribution by the Fixman term, which must be corrected for to 
unbiasedly sample the phase space. Our results show that this correction is warranted (Fig. Be) and (f)). 
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Both phase space sampling methods we introduce in this work are based on theoretically correct interpretations 
of the constrained partition function. The Monte Carlo phase space sampling method has the advantage of not 
requiring the Fixman correction to produce the correct average force prediciton. Additionally, this method may prove 
more useful for geometries such as polymer chains with a constant end-to-end distance. In such cases, selecting 
conformations that fit a given distribution is difficult but perturbing a conformation that already lies within that 
distribution is comparatively easy. The chief disadvantage of this method is the computational time required to 
perturb the chain, and the large number of steps required to adequately explore the conformational space. 

In contrast, the equilibrium phase space sampling, where uncorrelated conformations are selected, and momenta 
assigned after, is much quicker for free chains. This lack of correlation may more easily lend this method to large-scale 
parallelization as well. One minor disadvantage of this method is that it introduces a discrepancy in the mean force, 
unless the sampled conformations are weighted by the Fixman correction. 


V. CONCLUSIONS 

We introduced two different computational methods to sample microstates of a freely jointed chain in an unbiased 
manner. Unlike mainstream Monte Carlo methods that explore the conformational space only, our methods explore the 
complete phase space to obtain full distribution of any velocity-dependent microscopic variable in thermal equilibrium. 
We applied these methods to a terminally pinned chain in full space and half space to calculate constraint forces. The 
force distribution was non-Gaussian and asymmetric with both tensile and compressive forces. Most notably, the most 
likely force was smaller in magnitude than the mean force, which is more commonly considered both theoretically 
and experimentally. The constraint force exhibited little chain-length dependence beyond ~ 10 monomer lengths. 
The presence of a plane leads to a larger mean constraint force. This work can be extended to investigate fluctuating 
forces using more realistic polymer models and confinement. 
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(a) (c) 




FIG. 1. The origin of force by a constrained chain. A freely jointed chain (blue balls and sticks) is terminally pinned to a 
point (crosshair) and can move in full space or half space (dashed horizontal line), (a) Constraint (centrifugal) force. When 
freely jointed mass points of the chain move with some velocities (red arrows), the pin point will experience a force. The mass 
points can also exert force on the pin point when they recoil from the wall, (b) Thermodynamic (entropic) force. If the chain 
is pinned through a longer tether, more conformations become available. For example, the hrst mass point can explore a larger 
surface (compare blue and gray circles). Some conformations prohibited by the wall boundary can become acceptable as well 
when the tether is lengthened. This increase in conformational space leads to an increase in entropy, and thus a net force away 
from the surface, (c) Generalized coordinates in the global frame. The first three mass points of the chain are shown as black 
balls. The hrst three Cartesian coordinates Z\) dehne the position of the hrst mass point with respect to the origin. 

The tether length is ri. The positions of the rest are dehned by global zenith {9) and azimuthal {(j)) angles. 
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FIG. 2. Chain pinned in full space, (a) Distribution of kinetic energy at different points along the chain. The first mass point 
after the pinning position has less energy than average, and the last one has slightly more. The mean kinetic energy remains 
NksT in all cases. Shown are results from different chain lengths (A^=10, 20, and 30). (b) Distribution of radial constraint 
forces. While the mean constraint force equals that due to a single mass point, additional links in the polymer chain increase 
likelihood of strong pulling forces and even introduce the possibility of compressive forces. The most probable force is near 
— 1 shown by the vertical dotted line. Distributions are computed for different chain lengths {N=2, 9, and 29) from 150 K 
samples, (c) Mean constraint force and entropic force. The average force (black circles) is independent of length, and matches 
the predicted entropic force (red line), which results from the hrst mass point confined to the surface of a sphere moving to a 
larger radius, (d) Standard deviation of constraint forces. The standard deviation is on the same order of magnitude as the 
mean force, and quickly saturates as chain length increases. 
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FIG. 3. Origin of compressive forces. Some fraction of conformations (blue) and velocities (red) produce compressive forces 
along the first link in the chain, in contrast to the expectation from a single link acting as a pendulum. These occur when the 
chain is sharply bent, and the downstream mass points are moving more rapidly than those closer to the pin point {v 2 > 111 ). 
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f{xa(3) Length (/a) 


Length (/a) 


FIG. 4. Chain pinned to plane, (a) Distribution of forces for phase-space sampled chains. The surface effect shifts the mean of 
the distribution slighty, but the shape is not changed relative to the free case. The most probable force is near —1 shown by the 
vertical dotted line. Distributions are computed for different chain lengths {N=2, 9, and 29), from 150 K samples, (b) Mean 
force from chains sampled via phase space Monte Carlo (circles) and thermodynamic force from conformation space partition 
function (dashed line). The average force lies above the 2/a/3 value of the free chain case, as additional entropy is gained by 
pulling away from the surface. The mean force exhbits some length dependence not seen in the free chain case, however it 
quickly saturates around 10 links in the polymer chain, (c) Standard deviation of sampled forces. As in the free chain case, 
the deviation is on the same order of magnitude as the mean force. It grows slightly as a function of length, quickly saturating 
around 5-7 links, (d) Distribution of forces for conformational space sampled chains. The same distribution observed in (a) 
persists as we consider longer chains {N=2, 9, and 99). The most probable force is near —1 shown by the vertical dotted line. 
The curves in the figure represent 300 K samples, (e) Mean force from equilibrium sampled chains, before applying the Fixman 
correction (x) and after (o), displayed along side the partition function prediction (dashed line). The mean force remains 
constant as we explore larger lengths than were accessible in (b). (f) Standard deviation for force from equilibrium sampled 
chains, with (o) and without (x) the Fixman correction. 
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SUPPLEMENTAL INFORMATION 
A. Dependence of the mass metric on the end constraint 

The Cartesian coordinates of the n-th mass point are given by (n > 2). 

N 

Xn=Xi + E sin 6i cos (l)i 

i^2 

N 

Vn = yi + sin 9i sin (/)* 

i^2 

N 

Zn = Zi + cos Oi . 
i=2 

It is assumed that all adjacent mass points are joined by an inextensible linker of unit length. The spherical coordinates 
of the first mass point are given by 


Xi = ri sin di cos (pi 
jji = ri sin 01 sin pi 
zi = ri cos 6 * 1 . 


We define two sets of generalized coordinates u = {ri,6i,pi} and a subset of unconstrained coordinates q = {9i,pi} 
where i = 1,2,..., N. The Jacobian matrix J with respect to u can be constructed from smaller Jacobians for individual 
mass points. The Jacobian matrix is a 3 x 3 matrix for the first mass point: 


'dxi 

dxi 

dxi 

dri 

d9i 

dpi 

dyi 

dyi 

dyi 

dri 

d9i 

dpi 

dzi 

dzi 

dzi 

-dri 

Wi 

dpi- 


sin 9i cos pi ri cos 0i cos pi —ri sin 6 *i sin pi 
sin 6*1 sin ri cos 6 *i sin (/)i ri sin 6 *i cos (/)i , 

cos 01 —risin 0 i 0 


and a 3 X 2 matrix for the rest: 


dxi 

dxr 

d9. 

dpi 

dy^ 

dyi 

d9. 

dpi 

dzi 

dzi 

ld9. 

dpi- 


cos 9i cos pi 
cos 9i sin pi 
sin 9i 


— sin 9i sin pi 
sin 9i cos pi 


0 


The total Jacobian can be expressed as a block matrix: 


Ji 0 ••• 0 

Ji J 2 • • • 0 


Ji J 2 • • • 3n 


The mass matric becomes a {2N + 1) x (27V + 1) square matrix: 


TVJTJi (7V-1)JTJ2 


i-'i 

(7V-l)JjJi (7V-1)JJJ2 


'AfJl 


J^J2 


JJJat 




M„ = mJ^J = m 
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where ri-dependence is completely isolated in the first block row and column. Expanding the block matrices, one can 
obtain ri-dependence of the matrix elements: 


M„, = 


'N 

0 

0 

0 

Nr\ 

0 

0 

0 

Nrl sin^ S 


oc ri 

OC ri 


oc ri 

oc ri 


oc ri 
1 oc ri 


• • oc ri 

• • oc ri 


where oc ri denotes terms proportional to ri. We can obtain a 2N x 2N submatrix by excluding the first row 
and column of M„: 


Mg = m 


Nrl 


0 


OC ri 


0 Nrf sin^ 6*i oc ri • • • oc 


oc ri 
ri 


oc ri oc ri 


oc ri oc ri 


(SI) 


It is straightforward to show det(Mg) oc rf, which gives rise to an entropic force 2kBT/r upon constraining ri = r. 
To demonstrate this, we transform the partition function for the unconstrained system to the constrained one [30] . 
We use the Hamiltonian of a freely jointed chain without the end constraint to express the partition function: 

is a vector of 2N + 1 momenta conjugate to u. Delta functions are introduced to apply the end constraint on rii 

Pui 


Q,{r)= / 


where r-dependence of the constrained partition function is explicitly shown. The delta functions can be handled by 
decomposing the mass metric and using the transformation rule = d'^^Pq{Nmdri) 

Q^{r) = J - r)S{h)dri{Nmdfi)d‘^^qd‘^^Pq 

= NmJ 

Integrating with Pg, we can obtain the r-dependence of the partition function: 

Qdr) oc J det{Mq{r))d‘^^q (S2) 


B. Calculation of the metric determinant 


Chain conformations acquired from equilibirum sampling of the conformational space, without consideration of 
momentum space, will require a weighting factor given by the determinant of the metric tensor in the constrained 
space. The metric in the constrained space is generally dense, and so calculating the determinant g°‘ is computationally 
intensive. Fixman [18j made use of the simple form of the metric determinant g in the unconstrained space, to find 
the desired 5 “ from the determinant g^ in the smaller, orthogonal subspace spanned by the constrained dimensions. 

9 = 5 “/ ^ 5 “ = 5 // (S3) 

This is further simplified by considering the inverse H of the metric for the orthogonal subspace of constrained 
dimensions. This has the form 


N 


fc=0 


dri dvj 

dxk dxk 


(S4) 








15 


where is the bond length between points Xi_i and x^. The terms in the sum are only non-zero when i = j ± 1 and 
fc = i or i — 1. This ultimately produces a tridiagonal matrix 


■ 2 

-f 1 • f2 

0 

0 

0 

-f 1 • f2 

2 

-f2 - fa • • • 

0 

0 

0 

-f2 • h 

2 

0 

0 

0 

0 

0 

2 

-fAT-l • fw 

0 

0 

0 

-fiV-l • fw 

2 


Expanding this matrix in minors gives a recursive formula for the determinant h, where f„ is the unit relative position 
vector pointing from (n-l)-th mass point to n-th mass point. 

hji = (fn • ryj_i) hn—2 

For equal mass points, we begin with the base cases ho = hi = 1. Note this is slightly different than the Fixman case 
{hi = 2) as we are also constraining the position of the first mass point. Once we have computed h, the determinant 
is easily found from 


9^ = g-h 


(S7) 





